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We explore the dynamics of non-interacting particles loaded into a phase-modulated one- 
dimensional lattice formed by laterally oscillating square barriers. Tuning the parameters of the 
driven unit cell of the lattice selected parts of the classical phase space can be manipulated in a 
controllable manner. We find superdiffusion in position space for all parameters regimes. A directed 
current of an ensemble of particles can be created through locally breaking the spatiotemporal sym- 
metries of the time-driven potential. Magnitude and direction of the current are tunable. Several 
mechanisms for transient localization and trapping of particles in different wells of the driven unit 
cell are presented and analyzed. 
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I. INTRODUCTION 

Time-driven systems represent a major focus in sev- 
eral versatile research fields, such as the physics of atoms, 
molecules or mesoscopic systems [T}^. Generally, they 
are evoked by the occurrence of time-periodic forces. As 
an example, the dipole interaction of atoms exposed to 
laser fields gives rise to several interesting phenomena like 
laser stabilization, high harmonic generation or above 
barrier ionization. Another example is the coherent con- 
trol of quantum molecular dynamics by means of shaped 
femtosecond laser pulses. 

Among the most prominent mesoscopic systems are 
the various driven lattice setups, i.e. particles in an one- 
dimensional static potential are acted upon additionally 
by external time-dependent forces of zero mean. Exper- 
imentally they have been realized in condensed matter 
[TUl [TT] and cold atomic systems (see [T^HTl] and Refs. 
therein) . A remarkable observation is that these systems 
can show directed transport for an ensemble of parti- 
cles, although there exists no net force. Therefore, they 
are called ratchets. Originally the generation of ratchet 
effects has been addressed by employing external noise 
[3HS]- Similarly, the role of dissipation has been stud- 
ied thoroughly. In Ref. [7] it has been shown that di- 
rected transport occurs for underdamped particles in a 
sinusoidally rocked spatially asymmetric periodic poten- 
tial. Moreover, as the amplitude of the external driving 
is varied, the current flow is reversed several times. The 
underlying mechanism responsible for the existence of di- 
rected currents and the reversal of the transport has been 
identified in Refs. [3 [9]. Due to dissipation transport- 
ing attractors in phase space emerge and by choosing the 
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initial conditions appropriately it is possible to populate 
them selectively. The ratchet effect can also be generated 
in systems without dissipation and noise. In this case one 
speaks of "(deterministic) Hamiltonian ratchets". Re- 
cently, in Ref. [14 the directed transport of atoms in a 
BEG loaded into a flashing ratchet potential was demon- 
strated. The current flow appears only if certain tempo- 
ral and spatial symmetries of the driven potential, which 
have been identifled in Ref. [TS], are broken. The occur- 
rence of directed transport in Hamiltonian systems has 
been reported for the case of fully chaotic dynamics [TBl - 
[15] and for systems with mixed phase space [TOHIS]- For 
the latter case a sum rule for the transport veloc- 

ity of a classical ensemble of particles has been derived. 
Additionally, these authors have shown that in the semi- 
classical limit the quantum and the classical transport 
velocity coincide. As an extension in Ref. 126] the im- 
pact of avoided crossing between different transporting 
Floquet states has been considered. Tuning the control 
parameters leads to an enhancement or suppression of 
the current flow. Moreover, in Ref. [SJ the influence 
of an additional dc bias on the directed transport of a 
Hamiltonian ratchet has been studied. The authors have 
found the persistence of transporting invariant subman- 
ifolds like regular islands. In their vicinity trajectories 
can get sticky, such that they perform ballistic-like mo- 
tion. Remaining chaotic trajectories are accelerated by 
the bias field getting separated very fast from the ballistic 
type dynamics. 

For all previous setups the static potential is exposed 
to a so-called global driving law p!3lfT5l l2TI - [23] . i.e. the 
force acting upon the particles can be separated into two 
parts, which depend only on the spatial coordinate and 
the time, respectively. For the systems we explore in 
the present work this is not true anymore. Each po- 
tential barrier of the underlying lattice will be equipped 
with its own characteristic driving law, which gives rise to 
several new intriguing phenomena and the possibility to 
"locally engineer" the classical phase space. By adjusting 
carefully the parameters of the barriers and accordingly 



2 



their driving laws specific parts of the phase space can 
be manipulated in a controllable manner, whereas the 
remaining portion stays mainly unaffected. Moreover, 
the symmetries derived in 15^ can be broken by impos- 
ing spatially dependent phase shifts to the driving laws 
of the barriers, which will be called a "phase-modulated 
lattice" in the following. Thereby, a directed transport 
of an ensemble of particles is evoked and both the direc- 
tion and the magnitude of the current flow can be tuned 
easily. Importantly, for specific ranges of the barriers' 
potential height the particles show different localization 
behavior depending on their location within the unit cell 
of the lattice. 

The paper is organized as follows. In Sec. [ll] we give 
a detailed description of our model and define the setups 
which will be studied. In Sec. |III| the phase space is 
analyzed by means of stroboscopic Poincare surfaces of 
section (PSS). In this regard the desymmetrization of the 
phase space, meaning the loss of symmetry with respect 
to p ~ 0, is explored. Furthermore, we study the occur- 
rence of regular elliptic islands in the PSS, present an 
approximation for the velocity regime of the last stable 
torus and discuss the impact of cantori on the dynamics 
of trajectories. Section |IV| is devoted to the transport 
and localization properties of the setups and how they 
originate from the phase space properties of the consid- 
ered system. Finally, a conclusion and outlook is given 
in Sec. IVl 



II. SETUP 

The classical dynamics of an ensemble of identical, 
non-interacting particles in an one-dimensional, infinitely 
extended lattice of laterally oscillating square potential 
barriers of equal height Vq and width I is described by 
the Hamiltonian 



H{x,p,t) ^^ + V{x,t), 
where m is the mass of the particles and 



(1) 



\x-xo,, - Mt)\] (2) 



is the potential. XQ^i is the equilibrium position of the 
i-th barrier. fi{t) is a function with period T, which is 
explicitly allowed to depend on the site i. In the follow- 
ing, we will call fi{t) the driving law of the i-th barrier. 
The Hamiltonian ([ij is called to be unbiased, if the force 
F = —dV{x,t)/dx averaged over space and time van- 
ishes. This condition is fulfilled, when each barrier of 
the lattice is driven with a periodic driving law. Accord- 
ingly, a natural choice for fi{t) is a harmonic function. 
Therefore we set 
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FIG. 1: Schematic illustration of the lattice of laterally os- 
cillating barriers. 

where ipi G [0, 27r[ is a local phase shift, which depends 
on the i-th barrier. Here we have chosen the same am- 
plitude C and frequency uj for all barriers. Furthermore, 
we choose for the equilibrium position of the potentials 



Xo,: 



iD with > 0, 



(4) 



such that the static counterpart of Hamiltonian ([T]) 
is an equally spaced lattice of identical barriers, i.e. 
in their equilibrium position they are centered around 
{0, ±D, ±2D, . . .} (see Fig. [l]). Fig. [l] shows a sketch of 
the lattice including the relevant parameters. Only for 
very special choices of the phases, e.g. {(pi = 0} for all 
fi(t), neighboring barriers have equal spatial distance at 
a given snapshot. To avoid overlap of the barriers for ar- 
bitrary choices of {ipi}, the distance between the barriers 
must obey D>2C + l. 

The focus of this work are the transport properties of 
the system described by the Hamiltonian ([iJ. Let us 
therefore discuss under which conditions directed trans- 
port can occur. InRef. [15] . the authors have proven that 
the directed current vanishes if the equation of motions 
are invariant under certain transformations, i.e. if certain 
symmetries hold. These transformations change the sign 
of the velocity and as a result it is possible to construct 
for every trajectory a mirror image with the opposite sign 
for the velocity. In this case, averaging over a represen- 
tative set of trajectories yields zero mean velocity, i.e. a 
vanishing directed transport. Such transformations re- 
verse cither time or spatial coordinates together with a 
constant shift. Thus they take on the appearance 



Ta : X —X -\- Cx , t-^t + Ct, 

Tb : X X + Cx , t ^ —t + Ct- 



(5) 



The equations of motion belonging to the Hamiltonian 
([!]) are simply given by 



mx 



dV{x,t) 
dx ' 



(6) 



fi{t) = Ccos{ujt + (fii), 



(3) 



For the absence of directed transport for symmetry rea- 
sons it is therefore necessary that the potential V{x, t) is 
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invariant under one of the transformations ([s]). We re- 
mark that a vanishing current can also occur if the above 
symmetries do not hold, which is however a non-generic 
case and is encountered only for specifically chosen pa- 
rameter values. In the case of the uniformly oscillating 
lattice 

ifi — const. Vi (7) 
actually two such transformations can be identified. 



T : X ^ ~x , t — > t + TT, 
T' : X ^ X , t ^ -t. 



(8) 



Consequently, for the uniformly oscillating lattice no di- 
rected transport will occur. The same is true, if the local 
phase shifts of the barriers alternate, i.e. 



ipi if i is even 
(f2 else. 



In this case the transformation is given by 



T 



D , t 



(9) 



(10) 



In order to reduce the number of varying parameters, we 
set the frequency and the amplitude of the driving laws 
to w = 1 and C = 1. Moreover, the barriers' width and 
the equilibrium distance are fixed to I = 0.4 and D = 4.4, 
respectively. Finally, without loss of generality the mass 
of the particles is chosen m = 1. We therefore remain 
with two parameters: the barrier height Vq and the local 
phase shifts {<p,;}. Since our major interest are the lat- 
tice's transport properties, we focus on setups, for which 
the symmetries ([s]) are broken. This is implemented by 
spatially modulating the phase shifts {(pi} of the driving 
laws. In the following three different setups are studied 
in detail: 

(a) Firstly, the potential height is fixed to Vq = 0.16 
and the sites are equipped by a linearly increasing phase, 
i.e. a constant phase gradient is chosen. 



f^(t) = cos U + -27r 
' n 



(11) 



where the periods of the site-dependent phases (for rea- 
sons of brevity we call them "phase periods") are n = 
1, 3, 6, 10. For n = 3, the sequence of phase shifts {(fi} is 
{. . . , 0, ^, 0, ^, . . .}, whereas for n 1 a lattice of 
uniformly oscillating barriers is recovered. One can verify 
straightforwardly that it is not possible to define proper 
Ct and Cx, such that the potential is invariant under one 
of the transformations ([5| for phase periods greater than 
or equal to three. Hence a nonzero mean velocity could 
be expected for n — 3,6, 10. 

(b) For the second class of setups we keep Vq = 0.16 
and the above phase gradient with period three setup is 
perturbed (perturbed phase gradient), i.e. the sequence 
of phase shifts becomes {. . . , 0, ^ ± a, 0, ^ ± a, . . .} 




FIG. 2: (Color online) Stroboscopic PSS of the lattice of uni- 
formly oscillating square potentials. For a better illustration 
we have added in (a) the potential energy to the kinetic en- 
ergy of the particles, in order to avoid the discontinuities in 
the "regular" PSS (see (b)). 



the equidistant phase gradient is 



such that for a 
recovered. 

(c) A specific phase period three 
{. . . , 0, Yj), 1^, 0, 1^, . . .} is chosen and the global 
potential height Vq is varied ( broken phase gradient) . 



III. ANALYSIS OF PHASE SPACE 

This section is devoted to the analysis of the phase 
space of the Hamiltonian ^ for the above-provided se- 
tups. General characteristics of the PSS, like the appear- 
ance of regular islands and the chaotic sea, are discussed. 
Especially, main emphasis is placed on the comparison 
of the phase space properties of the uniformly oscillating 
lattice n = 1 (global driving, n — 1) and the setups with 
spatial phase modulation (local driving, n > 3) . Further- 
more, the stability of the phase space structures against 
perturbation of the phase pattern is studied. 



A. Poincare surfaces of section: Uniformly 
oscillating lattices 

In this section the phase space in the case of global har- 
monic driving is analyzed. Firstly, the method, how the 
Poincare surfaces of section for this and all the upcom- 
ing setups are obtained, is described briefly. Afterwards, 
the phase space is discussed and a kinematic approach 
to determine the border of the chaotic sea is presented. 
Finally, the flux through cantori and the transit time of 
trajectories is determined. 
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Since the particles move at constant velocity between 
collisions with the barriers' edges, the Hamiltonian flow 
can be described by an implicit two-dimensional map 
(for details see PF). To visualize the dynamics, we 
make Poincare surface of sections (PSS) of the phase 
space. Due to the time-periodicity of the Hamiltonian 
H{x,p,t + T) = H{x,p,t) (T = 27r), an area-preserving 
surface of section is obtained by taking stroboscopic 
"snapshots" of the (x,p)-plane at times t = riT with 
n G N. Additionally, the Hamiltonian possesses the 
translation invariance H{x -\- L,p,t) = H(x,p,t), where 
L is the length of the system's unit cell. Generally, L is 
not the distance between the barriers at their equilibrium 
positions, but for our chosen setups L is a multiple of D, 
i.e. L = nD, where n is the phase period. Thereby, the 
space coordinate can be restricted to the length L of the 
Hamiltonian's unit cell. The desired Poincare surface of 
section is then provided by the set of points 

M = {{x{t + kT) mod L,p{t + kT))\k e N}. (12) 

Obviously, since the time between successive "snapshots" 
of the (a:,p)-plane is equal, the PSS is done for a specific 
phase of the driving. It is therefore not surprising that 
the appearance of the surface of section is not universal 
but depends on exactly this phase. An example of how 
this manifests itself in the PSS is given after the next 
paragraph. 

In Fig. [2] (b) the PSS of the lattice with a phase pe- 
riod one, i.e. a lattice of uniformly oscillating barriers 
(n = 1,L=D) is shown. This Poincare surface of sec- 
tion possesses discontinuities (dashed box in Fig. [2](b)), 
which are due to the fact that the particles are either 
inside or outside the barrier when the snapshot of the 
(x,p)-plane is taken. The particles outside have Epot — 
and accordingly Epot = Vo inside. Since the potential is 
not smooth, this provides discontinuities in the PSS. For 
illustrative reasons it is therefore useful to add for those 
particles, which are at the moment of the snapshot inside 
the barrier, the potential energy Epot — Vq to the kinetic 
energy. This has been done in Fig. [2] (a) . The invariant 
curves are then continuous, yet the barrier can be seen as 
a blank squared region in the PSS. Nevertheless, to avoid 
the discontinuities all the subsequent PSS are presented 
in this way. 

Another eye-catching feature of the PSS is reflection 
symmetry with respect to p = 0. However, this is not 
universally valid, but it occurs only for specific phases of 
the snapshot of the (a;,p)-plane. As Fig. [2] (a) shows, the 
PSS is done at the moment the barrier arrives at one of 
its turning points, i.e. the i-th barrier is centered around 
xo^i = iD + C at this time instant. For the first unit 
cell i = 0, this yields a;o,o = 1 (blank squared region in 
Fig. [2] (a)). If the PSS is made for instance when the 
barrier has its maximum velocity, then this symmetry 
is absent. Yet, the dynamics of the lattice of uniformly 
oscillating barrier can be completely classified by means 
of the stroboscopic Poincare surface of section [23] • 

For large momenta \p\ ^ 1 the potential is negligible 




FIG. 3: (Color online) (a) Central periodic orbit {x « 
2.1, p ~ 1.4) for the resonance with winding number w = 2 in 
the lattice of uniformly oscillating barriers, (b) Periodic two 
bounce orbit between the first and second barrier in the case 
of the specific phase period three {...,0,-^,^,0,-^,^,...}. 

due to its finite height. In this limit, the intcgrable dy- 
namics of a free particle is recovered and the phase space 
is foliated by invariant curves, which are topologically 
equivalent to a torus. In the PSS these curves appear as 
straight lines, which stretch out over the whole unit cell. 
This simple part of the PSS is not displayed in Fig. [2] (a) 
and occurs for \p\ > 3. With decreasing momentum, the 
integrability is lost, but large regular domains remain. 
Many orbits still lie on deformed tori, which occur as 
curved lines in the PSS (\p\ « 2 in Fig. [2] (a)). Parti- 
cles on these invariant curves travel through the lattice 
in the direction of their initial momentum. Indeed the 
KAM theorem predicts that tori with sufficiently irra- 
tional winding number survive under a small perturba- 
tion. The winding number is defined as the limit 

^^^.^.(O-.(O)^ (13) 

if it exists, x and t are measured in multiples of the 
spatial period of the unit cell and time period, respec- 
tively. Thus w is proportional to the average velocity 
of the particle in the lattice w ^ v. Tori with rational 
winding numbers w = r/s are excluded in the KAM the- 
orem. According to the Poincare-Birkhoff theorem they 
dissolve into an even number of alternately elliptic and 
hyperbolic fixed points of period s. The trajectories in 
the extended system, which correspond to these periodic 
orbits in the PSS, travel r spatial unit cells in s time 
periods in the direction of their initial velocity. Hence 
in the PSS they occur as s distinct points. Around each 
elliptic periodic orbit there is again a set of invariant 
curves, which can be seen in our PSS as elliptic islands. 
Completely analogous to the trajectory belonging to the 
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periodic orbit in the center of this structure, the mo- 
tion of particles, which are inside these islands proceeds 
through the lattice only in one direction. Accordingly, 
the trajectories intersect sequentially the PSS at differ- 
ent islands of the corresponding chain of islands. In the 
vicinity of the hyperbolic fixed points there is an infinite 
number of homo- and heteroclinic intersections of the sta- 
ble and unstable manifolds, which yields horseshoe type 
dynamics and hence the presence of chaos. For large ki- 
netic energies these chaotic layers are too small to be vis- 
ible in the PSS. The magnification of the PSS (see Fig. 
|4]) for 1.75 < p < 1.95 shows such a separatrix region 
around p « 1.93. However, with decreasing momentum 
the strength of the perturbation and thereby the size of 
these layers increases until finally a large chaotic sea in 
the PSS for small momenta develops. Trajectories, be- 
longing to this chaotic sea in the PSS, wander diffusively 
through the lattice. 

At p = 0, there is an elliptic island, which has been 
discussed in detail in [5S] . It corresponds to trapped mo- 
tion in the scattering region of a single barrier. Moreover, 
this island exists only for a certain regime of the parame- 
ters (Vq, I). Due to the point-like interaction between the 
barrier and the particles it is possible to determine ana- 
lytically [5S] the position of these resonances in the PSS 
by exploiting the symmetry properties of the correspond- 
ing central periodic orbit. Other dominant elliptic islands 
in the PSS are the ones with winding number w — r, i.e. 
the central periodic orbits are trajectories, which travel 
r lattice sites during one oscillation period of the lat- 
tice. Fig. [s] (a) shows the trajectory, which corresponds 
to the central periodic orbit of the r — 2 resonance at 
X ~ 2.1, p « 1.4 in Fig. [5] (a), together with the two bar- 
riers in position space. Obviously the phases belonging 
to the four collisions are symmetric with respect to ^, 
(^4 = ~ ?1j ^3 = ~ such that the velocity V4 after 
transmission through the two barriers equals the initial 
velocity vq. In the appendix [A] a method to determine 
the positions of periodic orbits with few collisions in the 
PSS is described using the example of this w — 2-orbit. 

Let us now discuss the border regions of the PSS 
(1.75 < p < 2.0 in Fig. [2] (a)), where the large chaotic 
sea is bounded by tori, which have not been destroyed 
by the perturbation, i.e. the driving. In this regime 
of intermediate kinetic energy more chains of elliptic is- 
lands are visible in the PSS. For even larger momenta 
the chaotic sea is bounded by an invariant curve, which 
is called the "first invariant spanning curve" (FISC). Fig. 
|4] shows a magnification of this region of the PSS for pos- 
itive momenta with the FISC shown as a dashed line. 
For negative momenta there exists of course a FISC, too 
(Fig. [2| (a)). Ahhough the position of the FISC in the 
PSS is very hard to be determined analytically, one can 
derive at least an approximation to its position in veloc- 
ity space. We recall that particles on the FISC are fast 
enough to be transmitted, even if a collision takes place 
when the barrier moves with its extremal velocity u^ax 
in the same direction, i.e. the minimal relative kinetic 
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FIG. 4: (Color online) Magnification of the Poincare surface 
of section in the region of the upper boundary of the chaotic 
sea. The dashed curve is the last stable KAM torus, i.e. this 
invariant curve delimits the chaotic sea. 



energy must exceed the potential height. Thus we get 



(14) 



where v is the velocity of the particle. In case of a har- 
monic driving law, the extremal velocities of the barrier 
are given by u^x = ±ujC. Plugging u^x in equation (14) 



yields an approximation to the velocity regime where the 
FISC is located in the PSS. Depending on the sign of v 
we find 



m 



(15) 



where the index ± denotes the position of the FISC in the 
PSS for positive and negative velocities, respectively. In 
case of the harmonic driving law we have jf+j = in- 
dependent of the phase period. For m — 1,uj = 1, C=l 
and Vb — 0.16, we find from equation (151 \v±\ ~ 1.57. 
Fig. |4] shows that the FISC of the uniformly oscillating 
lattice is located indeed at p ~ 1.9 > m ■ v^. This dis- 
crepancy is significant (around 15 %), which is due to 
the fact that this approximation accounts only for a sim- 
ple kinematic consideration of the dynamics between the 
particle and a single barrier and neglects all the dynam- 
ical processes happening in the extended system. Nev- 
ertheless, this naive approach provides some insights, as 
we shall argue briefly in the following. For driving laws 
with more than one frequency, which have been used for 
instance in Ref. |23j . the extremal velocity u^x of the 
barrier depends on its sign. For example, in the case of 
the biharmonic driving law f{t) = C{sm{ujt) + cos(2wt) 
one finds u^x = 2ujC and Uex = — 3wC, respectively. In 
this case equation ( 15 ) yields v+ = 2u}C + \J 2 Vg / m and 
w_ — ~3u}C — \j2Vi^lm. Consequently, we expect that 
the FISC for positive momenta is located closer to p = 
than for negative momenta independently of the phase of 
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FIG. 5: (Color online) Poincare surface of section of a single 
chaotic trajectory, which has been started in the upper region 
and stopped once it has crossed the p = 0-axis. The diamonds 
and stars are periodic orbits corresponding to truncations of 
the continued fraction expansion of wi = and W2 = 

respectively. 



the PSS. Indeed in Ref. 
observed. 



1231 a similar behavior has been 



B. Flux through cantori and transit times 

Let us now turn to the discussion of the transit times 
of orbits in the PSS. First of all we give an example of 
trajectories in the chaotic sea, which are confined to a 
subpart of phase space for a very long time. Then a pro- 
cedure how to approximate the average escape time from 
this region is described briefly. Finally, this method is ap- 
plied exemplarily in the case of the uniformly oscillating 
lattice. 

In Fig. [5] a magnification of the PSS of a single trajec- 
tory is shown. It has been launched in the chaotic layer 
close to the FISC and once it has crossed the i-axis, the 
simulation has been stopped. Obviously, there are sud- 
den changes in the density of points in the PSS. This is a 
hallmark of the so-called cantori, which are remnants of 
dissolved tori with irrational winding number. These ob- 
jects can be regarded as tori with gaps, so that the Hamil- 
tonian flow is able to pass through [55]. The diamonds 
and stars in Fig. [5] are periodic orbits, which belong to 
truncations of the continued fraction expansions of can- 
tori with the noble winding numbers wi = ^^^^ — [2, 1°°] 

and W2 = = [2, 1, 1, 2, 1°°], where 7 = is the 

golden mean. The associated rational winding numbers 
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FIG. 6: (Color online) Flux through the chains of periodic 
orbits, belonging to convergents of w\ and u;2 , as a function of 
the level j. The black curve is ^r^/s^ ~ C^"^ with ^ ^ 4.339. 



for the two periodic orbits are 



V, = [2, 1^] - 



144 
~55 



., = [2,1,1,2,15] = 



129 
'50' 



(16) 



As long as the trajectories are confined in the PSS to the 
region above the cantorus, their velocity never changes 
sign. Moreover, the magnitude of the velocity varies only 
over a narrow interval. Thus, they perform ballistic-like 
motion during this time span. This has severe impact on 
the dynamics of an ensemble of particles, which will be 
discussed in the next section. In order to obtain the av- 
erage length of the ballistic flights, we have to determine 
the flux $^ through a cantorus with irrational winding 
number w. Mather |30) has shown that the sequence 
of differences in action of periodic orbits belonging to 
truncations of the continued fraction expansion of w con- 
verges to 



lini W, 



r I s 



T I S 



(17) 



where Wr/s, ^r/s action of the "minimizing" and 

the "minimax" orbit respectively [29] . i.e. periodic or- 
bits with winding number w = r/s, belonging to the 
minimum (saddle point) of the action. The minimiz- 
ing orbit is generically unstable and hyperbolic, whereas 
the minimax orbit is either an elliptic or hyperbolic- 
with-reflection periodic orbit. In the following '^r/s — 
Wr/s — '^r/s called as the flux through a chain of pe- 
riodic orbits with winding number w = r/s. Due to 
the point-like interaction, the action W^/s of a periodic 
orbit can be calculated very easily. Between successive 
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collisions with the barriers' edges, the particles move bal- 
listically in a constant potential. Hence the Lagrangian 
is simply 



L{x, i, t) = 



mv 
2 



particle between barriers 
Vq particle in barrier 



(18) 

Consequently, the action of a periodic orbit in the PSS 
of period s is given by 



E 

4=0 



L{x, i, t)dt 



(19) 



ti = i2-K is the moment, when the i-th PSS is taken. Since 
the interaction is point-like this reduces to 



S-l Jmax(») 



r/s 



E 

i=0 



r I s 



(20) 



with 




where At 



- t 



particle between bariers 
i,j particle in barrier, 

(21) 

are the discrete time inter- 



vals between the z-th and the 2+ 1-st intersection of the 
trajectory with the surfaces of section, i.e. ti^o = i27r, 
^*jmax(i) = (* + l)27r and tij(i) are the points in time of 
the collisions occurring between the particle and one of 
the barriers' edges during ti to i^+i. For the search of the 
periodic orbits a variational scheme has been adopted, 
which allows the detection of orbits with periods up to 
many thousands. In Appendix |B] a description of this 



method is presented. Fig. 



£1 



shows the flux 



Tjjs 



through convergents of w\ ana W2 as a function of the 
level 2 of the truncated continued fraction expansions. 
In the beginning the flux scales according to the power- 
law ^r jsi ~ G^-^ with i K, 4.339 [51] (black curve in 
Fig. [6|. Still the sequence converges rapidly with in- 
creasing j as Fig. [6] shows. Asymptotically, we find for 
the fluxes through the cantori — 3.62 • 10^^ and 
<I>^2 = 1.09 • 10~^. An approximation to the transit time 
for particles to get from the region above the cantorus 
with winding number Wi to the phase space below it is 
given by 



tit); 



(22) 



where A is the area in the PSS above Wi. By dividing 
the PSS into small squares, we have estimated these ar- 
eas to be = 0.12 and = 0.18, which yields for 
the total transit time ttransit = + = 3.48 • 10^ to 
get in the PSS from the region above the cantorus w\ 
to the region below the cantorus w^- Consequently, on 
average the particles get confined 348.000 periods of the 
driving in this region of the phase space and thus per- 
form during this time ballistic fiights. We note that the 



above transit time is meant to provide a rough estimate 
for the transport through the cantori. The comparison 
with numerical data obtained by simulating an ensem- 
ble placed in a chaotic region close to the FISC shows 
that the transit time is in general even longer. Beside 
the effect of being trapped to the region above the can- 
torus, the trajectories accessorily can get sticky to the 
hierarchy of elliptic islands surrounded by subislands be- 
ing above the cantorus. The latter process extends the 
ballistic flights signiflcantly. 



C. Phase-modulated lattices 

So far we have considered only the lattice of uniformly 
oscillating barriers (global driving), i.e. = const. Vi}. 
Let us now turn to the setups, which have been defined 
in section |ll] We will discuss especially the differences 
concerning their phase space structure compared to the 
setup with global driving. Firstly, the setup with con- 
stant phase gradient is studied. Then the stability of the 
phase space structure with respect to a perturbation of 
the gradient is considered. Finally, we study the impact 
of the potential height Vb on the dynamics. 

a. Constant phase gradient. In Fig. [t] (a)-(c) the 
PSS of the lattice with a linear, equidistant phase gra- 
dient of period 3, 5 and 10 is shown. Contrary to the 
case n — 1, the PSS is obviously not symmetric with 
respect to p = and moreover there is no alternative 
phase of a chosen PSS for which this symmetry is re- 
stored. As a result not only the PSS but the complete 
phase space of these phase periods is asymmetric with 
respect to p = 0. In the next section we will see that 
this desymmetrization is actually the origin of the occur- 
rence of directed currents in the system. Nevertheless, 
the PSS shows in general the same structure, which we 
have already discussed for the lattice of uniformly os- 
cillating barriers. For small momenta there is a large 
chaotic sea with embedded resonances. Since the poten- 
tial height Vo and the potential width / are equal for all 
phase periods, each barrier possesses analogous to n = 1 
at p = a small elliptic island of bounded motion in its 
scattering region [5S] (e.g. xi fa 2, X2 ~ 6, 2:3 w 11 for 
ri = 3 in Fig. [7] (a)). With increasing kinetic energy 
more chains of resonances appear until the chaotic sea is 
bounded by the FISCs. Ahhough the form of the FISC 
depends on the phase at which the PSS is taken, Fig. [7] 
shows that for all periods of the gradient it is located in 
the velocity range defined by equation (15 1. However, a 
significant symmetry the phase space possesses for n = 1 
has disappeared, namely the areas of the elliptic islands 
in the PSS with winding numbers w = r/s and w = —r/s 
are usually not equal anymore. For example in the case 
n — 3 the dominant resonance is the one with winding 
number w = 1/2 (p « 1, a; ~ 0.5 in Fig. [t] (a)). On the 
other hand there are three resonances with w = —1/2 
(p w -1, xi w 0.5, X2 w 3, cca ss 5 in Fig. [t] (a)). Ap- 
parently, their total area in the PSS is less than the area 
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FIG. 7: (Color online) Stroboscopic PSS of the lattice with a linear phase gradient of period 3 (a), 5 (b) and (10). The length 
of the unit cell increases with the period of the gradient. The dashed curves are the FISCs. 



of the single w = 1/2 resonance. Since these areas are 
preserved under the Hamiltonian flow, this is universally 
valid for all phases of the PSS. 

b. Perturbed phase gradient. As the next step, the 
equidistance of the phase gradient with period 3 is bro- 
ken up by imposing an additional phase shift a to the 
"middle" barrier. For a relatively large range of a one ob- 
serves that the PSS remains on a coarse scale unaffected 
by this perturbation in the sense that the position and 
the size of the significant resonances stays approximately 
the same. However, close to the FISC important changes 
occur. Therefore the upper part of the chaotic sea is 
magnified in Fig. [s] for a = (a) and a = -OMe. (b). 
Obviously, the former stable FISC in the lattice with the 
equidistant gradient dissolves, when decreasing a from 
to — 0.037r/3. Thereby, the proportion of the chaotic 
layer is increased considerably. However, at the position 
where the FISC is for a = in Fig. [s] (a) a cantorus 
remains in the PSS for a — — 0.037r/3, which represents 
a strong barrier, such that the Hamiltonian flow passes 
through it very slowly. As a consequence it takes long 
until a trajectory, which has been started with small mo- 
mentum in the chaotic sea, samples the accessible region 
of phase space above the cantorus. Consequently, for a 
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FIG. 8: (Color online) Magnification of the Poincare surface 
of section in the upper region of the FISC for two different 
phase gradients, (a) is the equidistant gradient of period 3 
and in (b) the phase of the "middle" barrier has been shifted 
by a = — 0.037r/3. The dashed curves are the FISCs. 
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finite time this cantorus acts in a similar way like the 
FISC, i.e. as it would be an impenetrable torus. Yet, in 
the limit of long simulation time the Hamiltonian flow is 
able to pass through. Therefore these cantori are usu- 
ally called "partial barriers" and possess a major impact 
on the transient dynamics and the long-term transport. 
Thus, one has to be careful, when estimating the trans- 
port velocity by simulating an ensemble of particles with 
small initial momenta in the chaotic sea. 

c. Varying the potential height Vq. Let us now 
consider the impact of a variation of the global 
potential height Vq on the phase space properties 
of a lattice with phase shift period three sequence 
{...,0, f^,ff,0, ^,...}. In Fig. [9] the PSS for four 
different values of Vq is shown. Ooviously, the island 
at p = being a property of scattering from a single 
barrier and corresponding to a localized dynamics has 
disappeared, because the parameters (Vo,0 are chosen 
such that its central periodic orbit has ceased to exist 
[55] . Still, for not too high values of the potential height 
Vq < 13.0 the chaotic sea is connected (see Fig. 9] (a)). 
Yet, with increasing Vq another class of localized dynam- 
ics arises in the system. As Fig. [9](b) shows, a separated 
chaotic sea between the first and the second barrier arises. 
Obviously, this part of phase space is confined by impen- 
etrable tori, i.e. the particles starting in this sea cannot 
escape from it and are thus trapped between the barriers. 
(Note that the "empty barrier" region in the PSS is not a 
confinement criterion like the invariant tori.) As the po- 
tential height is increased further, this localized chaotic 
dynamics appears between the other barriers, too. Em- 
bedded in these chaotic seas are elliptic islands, e.g. at 
X = 9.1, p = 3.2 in Fig. [9] (c). Intuitively, the occurrence 
of these localized chaotic seas is surprising, because one 
could think that with increasing potential height the par- 
ticles just have to acquire more collisions until they are 
finally fast enough to surpass Vq. Even more astonishing 
is the fact that apparently these regions in phase space 
do not appear simultaneously between all barriers but 
successively within increasing Vo- Nevertheless, this be- 
havior can be resolved straightforwardly. For small mo- 
menta p^/2 ^ Vq the particles exhibit exclusively reflec- 
tive collisions, i.e. they cannot penetrate into the bar- 
riers. In this regime the dynamics of the Hamiltonian 
^ is equivalent to the Fermi-Ulam-Model (FUM) of a 
particle bouncing between two infinitely heavy oscillat- 
ing walls. In the case of smooth driving laws the phase 
space of the FUM is not globally stochastic but possesses 
a FISC at momentum pi,, which prevents particles from 
gaining arbitrarily high momenta [32) . Consequently, lo- 
calized chaos in the case of the driven lattice will occur 
if the potential height is larger than the kinetic energy 
associated with this momentum Vq > p^/2. For mo- 
menta less than a certain value ps the phase space of the 
FUM is indeed completely chaotic. In the case of the 
previous setups Figs. [7) |8] with Vq — 0.16 and in Fig. 
[9] (a) the potential height Vq has been chosen smaller 
than the corresponding kinetic energy Vo < Ps/2- Con- 



sequently, all orbits of the large chaotic sea move diffu- 
sively through the lattice in these cases. For intermediate 
momenta \ps\ < \pr\ < \Pb\ the phase space of the FUM 
possesses a well known resonance structure, whose cen- 
tral periodic orbits correspond to trajectories similar to 
the one shown in Fig. [3] (b). When the global potential 
height is increased, these orbits and their surrounding 
elliptic islands appear in the PSS between two barriers 
prior to the blocking tori, since \pr\ < \pb\- Due to the 
different phase relations the barriers have with respect to 
each other, these confined chaotic seas do not occur si- 
multaneously between different pairs of barriers but one 
after another with increasing potential height (see Fig. 
|9](b)-(d)). Of course, there coexists a chaotic, diffusive 
dynamics of the particles in the lattice, too. 

Let us briefly summarize the different possibilities for 
trapping particles in one unit cell of the lattice. In the 
case of small potential heights, which have been used for 
instance in setup (a) and (b), this can be done by means 
of the elliptic island at p = corresponding to trapping 
in a single barrier. For larger values of Vq this struc- 
ture is gone, but islands belonging to trapped trajecto- 
ries between two barriers are present in phase space (see 
discussion in the previous paragraph). Furthermore, for 
Vo > 13 there is the possibility of trapping particles in the 
confined chaotic seas. Finally, we want to remark that for 
setups with very specifically chosen potential heights and 
phase gradients it is possible to achieve regular, localized 
dynamics, which extends over the scattering regions of 
several barriers. However, this is a case, which requires 
fine-tuning, since the associated elliptic islands are very 
tiny and for the parameter regime chosen in our setups 
this kind of dynamics does not occur. 



IV. TRANSPORT AND LOCALIZATION OF 
PARTICLES 



In this section the transport and the localization prop- 
erties of the chosen setups are studied. As we have seen 
in section |IT] directed transport can be ruled out for the 
case of a lattice of uniformly oscillating barriers. Yet, 
whether the other setups show directed transport can- 
not be judged prior to a more detailed analysis, because 
breaking the symmetries (see Eq. ([5| and Ref. [TS]) is 
necessary but not sufficient for the occurrence of directed 
transport. 



Transport properties in phase-modulated 
lattices 



We will deal with the question whether the setups (a)- 
(c) show a directed flow of particles. In the course of this 
study the stability of the current against perturbations 
of the gradient and the diffusion properties in the lattices 
are discussed, too. 
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FIG. 9: PSS for a lattice with a gradient of period 3 and four different global potential heights. The values are Vo — 2.0 (a), 
Vo = 13.0 (b), Vo = 22.0 (c) and Vo = 32.0 (d). 



a. Transport properties for constant phase gradients 
In the upper part of Fig. [lO] the evolution in time of the 
absolute value of the mean position for several 

periods of the phase pattern is shown. All particles have 
been started with small momenta in the chaotic sea close 
to the origin of the lattice, in order to avoid that some 
of them are initialized already in a small elliptic island 
belonging to a ballistic flight. Obviously, for n = 1 (lat- 
tice of uniformly oscillating barriers) no directed trans- 
port occurs, as we have expected because in this case the 
relevant symmetries are not broken (see Sec. |ll]). The 
ensemble average of the position fluctuates only around 
(x) = 0. In order to make this visible the data-points 
with {x){t) < {{x){t) > 0) for n = 1 have been plotted 
with stars (diamonds). A similar behavior is observed 
for {v){t) (lower part of Fig. 10). For n = 1 the mean 



velocity changes its sign several times correlated with a 
corresponding behavior of {x){t). Contrary, the other 
phase periods show clearly directed transport. After an 
initial transient i w 10^ — 10"* the mean position grows 
according to (x) (t) = Umcan ■ t, where Wmean is referred to 
as the transport velocity of the system. An asymptotic 
linear fit to the curves (x) (t) yields the mean transport 
velocities presented in the first row of Table |l] For n = 3 
and n — 6 these values for Wmcan are consistent with the 
behavior of the ensemble average of the velocity (v) (t) for 
large times (see lower part of Fig. 10). After the initial 
transient, {v){t) never changes sign again and fluctuates 
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Whether this is also true in the case of 



the phase period n = 10 cannot be judged ultimately. 



TABLE I: Transport velocities for the different phase peri- 
ods obtained by performing asymptotically a linear fit to the 
evolution of the ensemble average of the position (first row) 
and by averaging over the chaotic sea of phase space (second 
row). 



since the fluctuations of {v){t) are of the same order of 
magnitude as Wmoan- Nevertheless, Fig. [lO] shows clearly 
that after the transient t w 10'' the ensemble average of 
the position grows linearly with time. 

Let us now discuss the question of the origin for the 
appearance of directed transport. Due to the choice of 
the initial conditions, the ensemble is localized in the 
chaotic sea between the FISCs for every point in time. 
Thus, the convergence of the ensemble average of the 
velocity {v){t) to a non-zero value can be explained by 
means of an asymmetry of the chaotic sea with respect to 
p = 0, if it is assumed additionally that the dynamics is 
ergodic with a uniform invariant density |23j . Again, it is 
important to emphasize that this desymmetrization has 
to occur in phase space and not only for a certain PSS, 
which depends on the phase of driving at the moment 
of the snapshot. For Hamiltonian system with mixed 
phase space the proof that the dynamics is ergodic in the 
chaotic subsets is still an open problem [3 3) and has been 
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achieved only for very special systems [M]. Nevertheless, 
if we assume ergodicity of the dynamics in the chaotic 
sea, then the phase space average of the velocity vqs 
over the chaotic sea should coincide with Wmoan- vcs can 
be calculated by 



vcs 



1 f 

^cs J m 
cs 



dT, with dT = dx dpd^, 



(23) 



where f2c5 is the phase space volume of the chaotic 
sea and ^ is the phase of the driving. Due to the 
time-periodicity of the Hamiltonian H{x,p,t + 2tt) = 
H[x,p,t), the volume Vies can be written as Vtcs — 
2tt ■ Acs with Acs being the area of the chaotic sea in the 
PSS. Performing the integration over x and p yields hence 
the mean velocity of the chaotic sea in the PSS, which is 
determined for a fixed ^ by dividing the PSS into small 
rectangles and averaging over the cells that get visited 
by a single long chaotic trajectory. This yields a func- 
tion v-pssiO^ whose convergence is checked by enlarging 
the grid of the PSS. Finally, the mean velocity of the 
chaotic sea of the complete phase space is obtained by 
averaging over the phases of the PSS, i.e. 



vcs 



1 

2^ 



(24) 



Of course this scheme is equivalent to the procedure pre- 
sented in |23j . We remark, that for our system it turns 
out to be more efficient to evaluate numerically the inte- 
gral ( 23 ) than applying the sum formula of [23] , where all 



the winding numbers and areas of the significant elliptic 
resonances in the PSS have to be determined. Further- 
more, even for very long simulation times {t > 10^) a 
"drift" of a chaotic trajectory beyond an intact KAM- 
torus by accumulating numerical errors has never been 
observed in our system. For the lattice of uniformly oscil- 
lating barriers, our procedure gives vcs — 0, as expected. 
In the case of other phase periods, i.e. different from one 
Eq. (24) yields non-zero values for wc5, which are sum- 
marized in the second row of Table|Tl For n = 3 and n = 6 
the velocities obtained by performing asymptotically a 
linear fit to the mean position and averaging over the 
chaotic sea coincide very well. However, for n = 10 Eq. 



(p3|) predicts that the transport should be actually in the 

This 
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opposite direction than it is observed (see Fig. 
contradiction can be resolved straightforwardly. Accord- 
ing to our discussion of the PSS, there are cantori in the 
chaotic parts of phase space. Depending on the system's 
parameters like phase period or potential height, the flux 
across them can become arbitrarily small and therefore 
these cantori restrict trajectories to certain subparts of 
phase space even for long simulation times. In the pre- 
vious section their impact on the appearance of ballistic 
flights has been discussed. Yet, as we have indicated al- 
ready in the discussion of the setup with the perturbed 
phase period three, there is of course the opposite ef- 
fect, too. For trajectories starting with small momenta 



like our initial ensembles, the cantori prevent the parti- 
cles even for long times from sampling the (transporting) 
phase space beyond the cantori. Indeed the surfaces of 
section of phase period n — 10 ioi t < 10^ show that the 
region of the chaotic sea especially close to the FISC for 
positive momenta has not been visited by a single tra- 
jectory at all. Since this part of phase space corresponds 
to a positively valued drift, a transient transport in the 
opposite direction occurs in our case. In order to verify 
this explanation, a smaller initial ensemble has been sim- 
ulated for a longer simulation time. For this simulation a 
zero crossing of {x){t) is observed at ta ~ 1.3 • 10^. After 
the zero crossing the mean position {x){t) grows in the 



positive x-direction, which is predicted by equation ( 23 1 . 
Moreover, for t > tcr the region in the PSS close to the 
FISC for positive momenta gets sequentially visited by 
trajectories. 

According to the above, the origin for the occurrence 
of directed transport is simply the desymmetrization of 
the chaotic sea with respect to p = 0, i.e. for the 
phase periods larger than three the phase space volume 
of the chaotic sea with p > and accordingly p < 
are not equal anymore. An important manifestation of 
this asymmetry is that the areas of elliptic islands with 
w = r/s and w — —r/s are different in the PSS. Usually, 
this phase space asymmetry is achieved by applying a bi- 
harmonic driving law to the static system [T3Hl5l l2Tti23] , 
such that the potential V{x, t) of the Hamiltonian breaks 
both symmetries derived in |15| . However, from the dis- 
cussion of the FISC in the previous section we know that 
a kinematic approach to the single barrier dynamics ex- 
plains already the asymmetry of the chaotic sea of the 
phase space with respect to p — 0, because the FISC, 
depending on the direction of the propagation, is located 
in different ranges of the velocity. For the lattice with 
phase-modulated harmonically driven barriers this sim- 
ple consideration, independently of the phase period, is 
not enough to explain the asymmetry. 

b. Diffusion properties for constant phase gradients 



In Sec. HI it has been shown that every barrier possesses 
an elliptic island at p — for the parameter values (Vq, I) 
chosen for setup (a). Accordingly, the particles in the 
chaotic sea, which move diffusively through the lattice, 
can become sticky to this structure. They therefore re- 
main at the same spatial location for many periods of 
the driving. On the other hand the particles can per- 
form ballistic flights. These events originate either from 
phases of motion during which the particles are confined 
by partial barriers to regions of phase space with non- 
zero average velocity or they are the result of stickiness 
to elliptic islands corresponding to periodic orbits, which 
travel through the lattice in the direction of their ini- 
tial momentum. Usually, strong partial barriers are close 
to the FISCs, where transporting islands are found, too. 
Consequently, both effects reinforce each other in this 
part of phase space, i.e. trajectories which are confined 
by a partial barrier and get additionally sticky to elliptic 
islands contribute to the longest ballistic flights in the 
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FIG. 10: (Color online) In the upper part of the figure the 
absolute value of the mean position averaged over an ensemble 
of 10^ particles, which has been started with small momentum 
p < 0.1 in the chaotic sea as a function of time is shown for 
various phase periods. For n = & the evolution of the standard 
deviation has been plotted in the inset. Accordingly, the lower 
part of the Figure shows the mean velocity as a function of 
time. 




FIG. 11: Cumulative distribution function P{t) for the length 
of ballistic flights (phase period n = 6). Fitting a power law 
to the distribution yields an exponent of /i = 1.6. 



system. 

These trajectories can be detected numerically very 
easily. We look simply for time intervals with a minimum 
length, where the sign of the velocity does not change 
upon interaction between the particle and the barriers, 
i.e. we search for ballistic flights. The corresponding 



probability distribution of the lengths of ballistic flights 
obeys a power law p{t) ~ i"". However, in order to 
reduce statistical noise [3S], we do not plot a simple his- 
togram of the data, but calculate the cumulative proba- 
bility distribution function P(t), which is defined by 



P{t) 



p{t')dt'. 



(25) 



Fig. 11 shows P{t) in a double-logarithmic plot. Clearly, 
there is a power-law behavior of P{t) ^ over sev- 
eral decades. The exponent can be estimated by a least- 
square fit using a power law, which gives fi = 1.6. Thus, 
due to equation (25 1 the probability distribution p{t) 
follows a power law with the exponent ly = fi + I = 
2.6. From continuous-time-random-walk (CTRW) the- 
ory PXl [531 it is known that this interplay between 
ballistic flights and waiting times should yield anomalous 
diffusion in configuration space. In the inset of Fig. [lO] 
the variance a{x){t) for n = 6 is shown, which obviously 
grows for large times according to a{x){t) ~ . Perform- 
ing an asymptotic power-law fit yields for the exponent 
7 « 0.65. CTRW theory predicts that between 7 and fi 
the relation 2j — 3 ~ fi [371 ISH] holds. For our numer- 
ical data this is quite well fulfilled. A similar behavior 
is found for the other phase periods, too. In the asymp- 
totic time limit the variance follows a power-law and the 
exponent is always between 0.5 and 1, i.e. the system 
shows universally superdiffusion in configuration space. 

c. Transport properties for the perburbed / broken 
phase gradient Now we turn to the discussion of the 
transport properties of the lattice, where the phase pe- 



riod three with equidistant gradient is perturbed. Fig. 12 
shows the transport velocity of the system as a function 
of the phase shift a of the "middle" barrier. For these 
setups and all the upcoming ones the transport velocity 
has been calculated exclusively by evaluating numerically 
Eq. (231, because this is much more efficient in terms of 



computational time than simulating the long term dy- 
namics of a whole ensemble. The star in Fig. 12 marks 
the setup with a — 
PSS shown in Fig. 



— for which one obtains the 
pi (b). Using this example one can 
again demonstrate very vivid the impact of the tuning of 
the asymmetric phase space on the direction of the trans- 
port. The PSS of the setup with equidistant phase period 
three, i.e. a = 0, shown in Fig. [S] (a) corresponds to the 
system with transport velocity vcs{<^ = 0) = —0.0456. 
By perturbing this equidistance of the gradient the for- 
merly stable FISC dissolves, which makes parts of the 
phase space volume with larger positive momenta acces- 
sible to chaotic trajectories. At the same time, the re- 
maining phase space, in particular the part with p < 0, 
remains basically unaffected, i.e. overall the phase space 
volume of the chaotic sea with p > Q increases signifi- 
cantly. Therefore, the phase space average of the velocity 
over the chaotic sea in equation ( 23 1 and thus the direc- 
tion of the transport changes its sign. For some values 
of a the directed transport vanishes, i.e. wcs(a) = 0. 
Yet, this does not imply that a symmetry of type ([5]) is 
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Transport velocity Umcan , as a function of the phase 
— ... of the "middle" barrier. With a star we 
have marked the setup, which we have chosen in section [Till 
to produce the PSS shown in Fig. |8](b). 



present. In fact both symmetries are broken. Neverthe- 
less, the phase space volumes with p > and p < of 
the chaotic sea are accidentally equal for some values of 
a. 

Finally, we consider the transport properties of the 
setup with variable potential height and fixed phase gra- 
dient of period three {. . . , 0, ff , 0, ^, ff , . . .}. By 
evaluating Eq. ( 23 1 the transport velocities depending 
on the global potential heights are calculated and sum- 
marized in Table [nj Similarly to the previous setup, the 
directed transport can be tuned by varying a parameter, 
which is now the potential height. 



B. Localization and trapping 

In this section the three different possibilities for trap- 
ping, which have been summarized in Sec. |III C[ are 
discussed in more detail. Especially the impact of the 
potential height is addressed, because the localization of 
particles is most sensitive to this parameter. 

a. Localization properties for a broken phase gradient 
In the course of the discussion of the diffusion properties 
it has been shown that the trajectories, which contribute 
to the directed current in the lattice, can still obey for 
long time phases of motion during which they are local- 
ized in the scattering region of one barrier. For small 
potential heights these events originate from stickiness 
to the elliptic island at p = 0, whose properties are de- 
fined exclusively by (Vo,0- Since these parameters are 
equal for all barriers, i.e. a single barrier property |28j . 
the stickiness of trajectories to this structure is identi- 
cal independently of the phase periods: The particles 



TABLE II: Transport velocities for different global potential 
heights Vb obtained by averaging over the chaotic sea of phase 
space. 



get localized in the scattering region of the first, second 
etc. barrier on average for the same time. In the case 
of very large potential heights, e.g. Vq = 32 for setup 
(c), there are confined chaotic seas between the barriers. 
Yet, the particles belonging to the directed flow cannot 
enter these regions in phase space, because these parts 
are separated by impenetrable tori. In the following we 
will show that it is possible for certain parameter val- 
ues to obtain a dynamics that exhibits phases of trapped 
motion in certain wells, i.e. in between certain definite 
barriers, of the driven unit cell of the lattice. To this 
end we look in the case of setup (c) for time intervals 
during which the particles, that move diffusively through 



the lattice, are localized between two barriers. Fig. 13 
shows for Vq = 4 the cumulative probability distribution 
of dwell times P{t) for particles between two consecutive 
barriers in one spatial unit cell, i.e. i = — 1 to i = 0, 
« = to i = 1 and i = 1 to i = 2, which will be labeled 
as well one, two and three, respectively. For this value of 



the potential height (see Sec. Ill C ) the elliptic island at 
p — has disappeared and thus there is no stickiness to 
this regular structure anymore. Furthermore, no trapped 
chaotic dynamics between barriers is possible. Instead all 
chaotic trajectories in the large chaotic sea move diffu- 
sively through the lattice. Obviously, P{t) is not equal 
for the different wells. A particle, which has been local- 
ized in the first well is trapped with less probability (for 
a time span t > 10'^) compared to a particle trapped in 
the third well. By recalling the discussion of the PSS for 
this values of Vq in Sec. |III C| this behavior can be re- 
solved. Accordingly, with increasing potential height the 
particles between two barriers are able to probe parts of 
the phase space of the corresponding FUM, which are 
located at higher momentum. However, the phase space 
of the FUM in general is very sensitive to the phase re- 
lations the walls have with respect to each other. Not 
only the position in momentum space of the FISC pb, 
but also both the momentum ps below which the FUM is 
completely chaotic and the momentum of elliptic islands 
belonging to trajectories shown in Fig. [3](b) depend crit- 
ically on the phase relations between the barriers. Due to 
the non-equidistant gradient these phase relations are not 
equal, i.e. the barriers belonging to different wells pos- 
sess different relative phases with respect to each other. 
Thus, for a fixed potential height it is possible that in 
a certain well only the chaotic part of the corresponding 
FUM's phase space is accessible to the particles, whereas 
in another well already elliptic islands embedded in the 
chaotic sea are present. Escape rates from these differ- 
ent parts of phase space are understandably not equal, 
because stickiness to regular structure leads to phases of 
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FIG. 13: (Color online) Cumulative distribution function 
P{t) of the dwell time for the three different wells with fixed 
potential height Vb = 4 for setup (c). In the inset P{t) be- 
longing to the second well for different values of Vo = 4 ... 8 
is shown. 



motion during which particles are trapped for long times 
between two barriers. By changing the phase gradient 
for a fixed potential height, it is possible to "engineer 
the phase space of the lattice in such a way" that the 
events with very long dwell times take place only in spe- 
cific wells. Of course, trajectories can also stay for a 
comparatively long time in a well, which does not con- 
tain any elliptic island at all, which is due to the so- 
called low velocity peaks [39j . Particles, that are slightly 
faster than the barrier at a collision, get decelerated to 
a very small velocity but still leave the scattering re- 
gion without experiencing a second collision. Afterwards 
they travel the spatial distance between the scattering 
regions of two barriers, which yields for the dwell time 
td = ■ Yet, the portion of trajectories, which 

are "trapped" according to this mechanism, is negligible. 
For Vq — A all wells contain elliptic islands. In fact, the 
characteristic asymptotic power-law behavior P{t) ^ 
of the cumulative probability distributions of dwell times 
is a hallmark of the associated stickiness to this regular 
structure in phase space. Generally, whenever a region 
in phase space contains elliptic islands, the occurrence 
of long algebraic tails in the escape rate is generic ^29]. 
Indeed, we observe that the exponent of the probability 
distribution p{t) ^ i"", which is related to P{t) through 
equation (25), is always larger than two. Similarly to 



the previous setups, the dynamics shows the characteris- 
tic interplay between ballistic flights and waiting times. 
Accordingly, we find that there is superdiffusion in coor- 
dinate space even for very large potential heights. 

b. Impact of potential height on localization properties 
Finally, let us discuss the variation of Vq infiuences P{t) 



for a given well. In the inset of Fig. 13 the cumulative 
probability distribution of the dwell time in the second 
well for various values of Vq between 4 and 8 is shown. 
With increasing Vq this is the first well with a localized 
chaotic sea (see Fig. |9] (b)). Especially for small and 
intermediate dwell times (10^ < t < 3 ■ 10"^) the function 
P{t) changes significantly. P{t) decreases less rapidly 
for larger values of Vb, which is understandable, because 
with increasing potential height the particles need to ac- 
cumulate more collisions in order to become fast enough 
to surpass the potential. Still, by doubling Vq the proba- 
bility for the particles to be trapped in the well for dwell 
times longer than 10^ grows up five orders of magnitude. 
Furthermore, the PSSs of trajectories with dwell times in 
this range show that beside the particles, which are sticky 
to elliptic islands, a great portion still obeys chaotic dy- 
namics. In order to answer the question how chaotic 
particles, which move diffusively through the lattice, can 
be trapped between two barriers for such a long time, 
we have to exploit one more time the FUM. Close to the 
FISC of the FUM are again cantori with small flux across 
them. Depending on the properties of these partial barri- 
ers the particles obeying chaotic dynamics are prevented 
for a significant time from sampling the phase space of 
the FUM, which is located at high momenta close to the 
FISC and associated to the kinetic energy a trajectory 
has to gain in order to leave the well. Consequently, 
once the potential height is chosen such that this cantori 
are stabilized, the escape rate from the region between 
barriers is lowered significantly. 



V. CONCLUSIONS AND OUTLOOK 

We have explored the dynamics of non-interacting par- 
ticles in one-dimensional, phase-modulated driven lat- 
tices. Depending on the parameters like the potential 
height and phase period the phase space of the system 
has been analyzed in detail. The impact of cantori on the 
transient transport properties, as well as the occurrence 
of ballistic flights due to the dynamical conflnement of 
trajectories to volumes in phase space with non-zero ve- 
locity, has been studied. Depending on the phase period 
a directed current is observed, whose occurrence has been 
traced back to the desymmetrization of the chaotic sea 
P5M25j . Commonly, this asymmetry of the phase space 
is achieved by breaking two established spatiotemporal 
symmetries [15j with an appropriately chosen biharmonic 
driving [TBHT51 [2TH25] . For this driving law a kinematic 
argument considering the single barrier dynamics is suf- 
ficient to explain the occurrence of a desymmetrization 
of the extended system's phase space, i.e. each barrier is 
transporting itself and the lattice inherits this property. 
In this work it has been shown that a simple harmonic 
driving law together with a local symmetry breaking im- 
plemented by local phase shifts is sufficient to generate a 
ffow of particles although each individual barrier is non- 
transporting. Consequently, the occurrence of the parti- 



15 



cle current is in this case a collective phenomena of the 
complete lattice. The direction and the magnitude of the 
transport can be tuned by the parameters of the system 
like the phase period or the potential height. 

Of course the most simplest way to revert the direction 
of the particle current is by inverting the gradient. With 
increasing potential height elliptic islands belonging to 
a trapping dynamics between the barriers arise. Their 
properties like size and position are influenced by the 
phase relation of the barriers with respect to each other. 
Accordingly, the phase space can be locally manipulated 
in a controllable manner, such that each well possesses 
its own characteristic escape rate. Furthermore, it has 
been shown that the system shows universally superdif- 
fusion in coordinate space. In retrospect this occurrence 
of anomalous diffusion is understandable, because its pre- 
requisites are always present in phase space. Independent 
of the system's parameters like potential height, phase 
period etc. the phase space is mixed. Thus, it possesses 
the hierarchy of elliptic islands surrounded by subislands. 
Events of stickiness to this regular structure yield either 
ballistic flights or waiting times. In both cases we have 
found that the corresponding probability distributions 
follow asymptotically typical power-laws, which give rise 
to anomalous diffusion and can be modeled by continuous 
time random walk theory. A major advantage of the pre- 
sented setup compared to other driven lattice systems 
is that the driving can be locally adjusted to engineer 
the phase space. Thereby, the desired transport and lo- 
calization properties can be achieved much more simply 
than by applying just another driving law. Furthermore, 
the local driving laws offers the possibility of building up 
"blocks of sublattices" , such that a current of particles 
can be manipulated. 
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Appendix A: Determining the position of periodic 
orbits in the PSS with few collisions 



In the following we describe briefly our method to de- 
termine the position of the periodic orbits with winding 
number w — r in the PSS using the example of w = 2. 
For this trajectory the initial velocity vq and the parti- 
cle's velocity between the first and the second barrier V2 
(see Fig. [s] (a)) are defined by the condition that the 
phases at collision with the first and third barrier are 
the same. In the case of the uniformly oscillating lattice 
{ifi = const. Vi}, this yields after some algebra 



_ L-l + 2cos(Ci) 



V2 = 



^ + 26 
L - I - 2cos(6) 



-26 



(Al) 
(A2) 



Due to the symmetry of the phases f upon collision, the 
particle's velocity after the first and the third collision 
with one of the barriers' edges are equal, i.e. vi = v^. 
Furthermore, the spatial distance between the collisions 
at phase 6 ^-nd 6 has to be covered by the particle during 
At = ^2 — ^1- Both conditions yield a system of two 
nonlinear equations 



/i(6, 6) = sin(ei) - sin(6) - v'K(6)+sin(a))'-2yo + \/(«2(6) + sin(6))' - 2^0 = 0, 
/2(ei,6) = cos(ei) - cos(6) + (6 - Ci) (voi^i) + sm(6))' - 214, - sin(6)) = 0. 



(A3) 



Solving equation (|A3|) for ^ and 6, the position of the 
elliptic island in the PSS can be calculated straightfor- 
wardly. In general, one can proceed similarly for setups 
with larger phase periods and for every resonance. A 
system of equations will therefore result 



We remark that it is difficult to generalize this scheme to 
periodic orbits with many collisions, since the approach 
to set up the system of equations and the roots corre- 
sponding to physical solutions depend on the specific 
symmetry properties of the phases upon collisions of the 
central periodic orbit belonging to the resonance. 



- ,a)=0, i = l...k}. (A4) 
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Appendix B: Finding periodic orbits with arbitrary 
winding number 

For periodic orbits with many collisions, which are 
needed for the calculation of the flux through a cantorus, 
the scheme presented in appendix |A] is obviously not fea- 
sible anymore. In this case it makes more sense to make 
use of a variational method, which will be presented in 
the following. Since the particles move ballistically, the 
trajectories in the lattice belonging to the periodic orbits 
in the PSS are completely defined by the moments of the 
collisions tk with the barriers' edges. Thus the action 



i.e. the action between successive collisions is added 
piecewise. N is the total number of collisions. Now we 
restrict ourselves to periodic orbits with winding num- 
ber w — r/s. Firstly, this yields that in equation (Bl) 
t^+i — ti + 27rs. Secondly, r specifies how many spatial 
unit cells a trajectory travels in the lattice, at which the 
velocity never changes sign. Accordingly, the number of 
collisions N, the phase period n and r are related through 
N — 2 ■ n ■ r. Inserting the Lagrangian of the system ( 18 ) 
in equation (Bl) gives 



( 20 ) of such a periodic orbit can be rewritten as 



N 

r/s — 

k=l 



Wk{tk,tk+l), 



(Bl) 



Wk{tk,tk + l) = { „fJt*'=+r*^f) ^^2 " (B2) 

2{tt+i-tk) ^0 (tk+1 - tk) particle in barrier 



particle between barriers 



Of course x(tk) equals with the position of one of the 
barriers' edges at a collision, which depends on the driv- 
ing law fi{t) at the i-th site. For particles in a barrier 
the spatial distance between the successive collisions is 



x{tk+i - x{tk)) = I + fi{tk+i) - fi{tk) and accordingly 
x{tk+i - x{tk)) = D - 1 + f,+i {tk+i) - Mtk) between the 
barriers. Setting m = 1 the action becomes 



I (-P ^+A+i(*fc+i) f^{tk)) particle between barriers 
Wk{tk,tk+^) = {a^f^^^^^l^^,^ (B3) 

2(tk+i—tk) 



ii+fdtk+i)-Mtk)) _ ^^^^^ _ particle in barrier 



We search for trajectories for which the action is 
extremal. The action gradient vector is VVF = 
{dW/dh, dW/dt^f with 

dU dt, dU ' 

For finding roots of VM^ a multidimensional Newton 
scheme has been applied. Therefore, the derivative ma- 
trix ("Hessian") of the action gradient is needed. It is 
a cyclic, tridiagonal matrix of the second derivatives of 
Wi, whose rank increases with the period s of the orbit. 



To get an initial starting point we start from the inte- 
grable limit (Vq = 0) and trace the periodic orbits as the 
perturbation, i.e. Vq, is increased. Finally, the stability 
of the obtained periodic orbit can be determined via the 
eigenvalues of the Hessian. In the case of the minimiz- 
ing orbit it has only positive eigenvalues, whereas for the 
minimax orbit there is a single negative eigenvalue. For 
orbits with not too high periods (s < 300) the results 
of this scheme have been compared additionally to other 
globally convergent methods [101 HI] and we have found 
good agreement. 
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